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ABSTRACT 

Variable 7-ray emission has been discovered in five Radio-loud Narrow Line Seyfert 1 (NLSy 1) galaxies by the 
Large Area Telescope (LAT) onboard the Fermi Gamma-Ray Space Telescope. This has clearly demonstrated 
that these NLSyl galaxies do have relativistic jets similar to two other cases of 7-ray emitting Active Galactic 
Nuclei (AGN), namely blazars and radio galaxies. We present here our results on the multi-band analysis 
of two 7-ray emitting NLSyl galaxies namely PKS 1502+036 (z = 0.409) and PKS 2004-447 (z = 0.240) 
using archival data. We generate multi-band long term light curves of these sources, build their spectral energy 
distribution (SED) and model them using an one zone leptonic model. They resemble more to the SEDs of 
the flat spectrum radio quasar (FSRQ) class of AGN. We then compare the SEDs of these two sources with 
two other Fermi detected AGN along the traditional blazar sequence, namely the BL Lac Mrk 421 (z = 0.03) 
and the FSRQ 3C 454.3 (z = 0.86). The SEDs of both PKS 1502+036 and PKS 2004-447 are found to be 
intermediate to the SEDs of Mrk 421 and 3C 454.3. In the 7-ray spectral index v/s 7-ray luminosity plane, 
both these NLSyl galaxies occupy a distinct position, wherein, they have luminosity between Mrk 421 and 
3C 454.3, however steep 7-ray spectra similar to 3C 454.3. Their Compton dominance as well as their X-ray 
spectral slope also lie between Mrk 421 and 3C 454.3. We argue that the physical properties of both PKS 
1502+036 and PKS 2004—447 are in general similar to blazars and intermediate between FSRQs and BL Lac 
objects and these sources thus could fit into the traditional blazar sequence. 

Keywords: galaxies: active — gamma rays: galaxies — quasars: individual(PKS 1502+036, PKS 2004-447 — 
galaxies: Seyfert 



1. INTRODUCTION 

The extra-galactic sky in the high energy 7-ray band 
is dominated by a special class of Active Galactic Nuclei 
(AGN) called blazars (Abdoetal. 2010b). These blazars 
which come under two sub-classes namely Flat Spec- 
trum Radio Quasars (FSRQs) and BL Lacs (see, for e.g., 
Urry & Padovani 1995), have similar broad band contin- 
uum properties, such that they emit variable non-thermal 
emission over the entire electromagnetic spectrum. They 
both have flat radio spectra (a r < 0.5; S v oc v~ a ), at 
GHz frequencies, high optical polarization, exhibits super- 
luminal pattern at radio-frequencies and also show rapid 
flux and polarization variations (Wagner & Witzel 1995; 
Andruchow et al. 2005). However, the distinction between 
FSRQs and BL Lacs is made on the difference in the strength 
of their optical emission lines. BL Lacs have weak optical 
emission lines with rest frame equivalent widths less than 
5 A(Stocke et al. 1991; Stickel et al. 1991). The observed 
broad band emission of blazars is believed to originate from 
relativistic jets which are aligned closely to the observer 
(Begelman et al. 1984). Their spectral energy distribution 
(SED) shows two prominent and distinct peaks: the low 
energy peak due to synchrotron emission, lies between the 
far-infrared and the X-ray band, and the high energy peak 
due to inverse Compton (IC) radiation, lies in the MeV to 
GeV energy range. The seed photons for IC emission can be 
synchrotron photons itself (Synchrotron Self Compton, SSC; 
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Konigl 1981; Marscher & Gear 1985; Ghisellini & Maraschi 
1989) or photons external to the jet (external Compton, 
EC; Begelman & Sikora 1987; Melia & Konigl 1989; 
Dermeretal. 1992). The plausible sources of seed pho- 
tons for EC can be accretion disc (Dermer & Schlickeiser 
1993; Boettcher et al. 1997), broad line region (BLR) 
(Sikora et al. 1994; Ghisellini & Madau 1996), dusty torus 
(Blazejowski et al. 2000; Ghisellini & Tavecchio 2008), 
bulge of the host galaxy, or cosmic microwave background 
radiation (Ghisellini & Tavecchio 2009). In the 7-ray regime, 
FSRQs are more luminous and exhibit softer spectral index 
than BL Lacs (Abdo et al. 2009a). It has been noted by 
Fossatietal. (1998) and Ghisellini et al. (1998) that both 
FSRQs and BL Lacs follow the so called blazar sequence. 
However, recently Giommi et al. (2012) has argued that the 
existence of blazar sequence could be due to selection effects. 

In addition to blazars which are dominating the 7- 
ray extra-galactic sky, Fermi has also discovered vari- 
able 7-ray emission from five Narrow Line Seyfert 1 
(NLSyl) galaxies (Abdoetal. 2009b,c; D'Ammando et al. 
2012; Calderone et al. 201 1). A few more NLSyl sources are 
suspected to be 7-ray emitters, although the detection signifi- 
cance is still low (Foschini 201 1). NLSyl galaxies are a sepa- 
rate class of AGN with peculiar properties. Their optical spec- 
tra are similar to the conventional broad line Seyfert 1 galax- 
ies except that they have narrow Balmer lines (FWHM (H^) 
< 2000 km/sec), weak [O III] ([O Ill]/Ha < 3) and strong 
optical Fe II lines (Osterbrock & Pogge 1985; Goodrich 
1989). They also have steep soft X-ray spectra (Boiler et al. 
1996; Wangetal. 1996; Leighly 1999a) and show rapid 
X-ray flux variations (Pounds et al. 1995; Leighly 1999b). 
These peculiar observational characteristics are attributed to 
them having low-mass black holes (~T0 6 — 10 8 M0) and ac- 
creting close to the Eddington limit (Peterson et al. 2000; 
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Hayashida 2000; Grupe & Mathur 2004; Zhou et al. 2006; 
Xu et al. 2012). However, recently, Calderone et al. (2012), 
using multi-wavelength data, have shown that radio-loud 
NLSyl (RL-NLSyl) galaxies have black hole masses simi- 
lar to blazars. These RL-NLSyl galaxies comprising about 
7% of NLSyl galaxy population, exhibit compact core-jet 
structure, fiat/inverted radio spectra, high brightness tempera- 
ture and superluminal motion (Komossa et al. 2006; Doi et al. 
2006). Recently, kiloparsec-scale radio structures are found 
in six RL-NLSyl galaxies (Doietal. 2012). The SEDs of 
some RL-NLSyl galaxies are found to be similar to those of 
high frequency peaked BL Lac objects and some of them also 
show a hard X-ray component in their bright optical/UV state 
(Yuan et al. 2008). Though these characteristics are indica- 
tive of the presence of relativistic jets in these sources that are 
closely aligned to the observer, the detection of 7-rays in five 
of them has confirmed their presence (Abdo et al. 2009b,c; 
D'Ammando et al. 2012) similar to the blazar class of AGN. 
Also, it has been recently reported that the optical and infra- 
red (IR) flux variations of some of these sources are simi- 
lar to blazars (Liu et al. 2010; Paliya et al. 2013; Jiang et al. 
2012; Maune et al. 2013). However, these similarities are at 
odds with the common belief that NLSyl galaxies are hosted 
in spiral galaxies with low mass black holes, while, blazars 
with powerful relativistic jets are hosted in elliptical galax- 
ies with high mass black holes. Therefore, detailed studies of 
these sources will enable us to understand the physical prop- 
erties of relativistic jets hosted in AGN with low mass central 
black holes and high accretion rates. In this work, we have 
carried out a detailed analysis of two such sources, namely, 
PKS 1502+036 and PKS 2004-447 to see (i) the similarities 
and/or differences of them vis-a-vis the blazars detected by 
Fermi and (ii) their place in the traditional blazar sequence. 
The detailed analysis of these sources are presented here for 
the first time. Throughout the work, we adopt __ TO =0.27, 
0/1=0.73 and Hubble constant H n = 70 km s _1 Mpc -1 . 

2. SAMPLE 

Among the 5 known 7-ray emitting NLSyl (7-NLSyl) 
galaxies, we have selected two of them, namely, PKS 
1502+036 and PKS 2004-447. They have intermediate 
redshifts when compared to other 7-NLSyl galaxies. PKS 
1502+036 has a very compact core on observations with the 
Very Large Array (VLA; Orienti et al. 2012; Kimball et al. 
2011). However, when observed at 15 GHz, using Very 
Long Baseline Array (VLBA) at parsec scale resolution, it 
was found to have a core-jet structure (Orienti et al. 2012). 
This NLSyl galaxy, with a central black hole mass of ~ 
4 xlO 6 M©, has a brightness temperature of ~10 12 K 
(Yuan et al. 2008). It is a faint, but persistent 7-ray emit- 
ter (Calderone et al. 201 1; Abdo et al. 2009c) and also show 
intra-night optical variability (INOV; Paliya et al. 2013) and 
rapid IR variability (Jiang et al. 2012). 

PKS 2004-447 has been classified as a NLSyl galaxy 
by Oshlacketal. (2001) and Gallo et al. (2006) but 
Komossa et al. (2006) classify it as narrow-line radio 
galaxy owing to the presence of weak Fe II lines in its 
optical spectra. It has a steep spectrum above 8.4 GHz 
(Orienti et al. 2012), however, has a flat spectrum below 5 
GHz (Healey et al. 2007). At 20 GHz, it shows a fractional 
polarization of 3.8% (Murphy et al. 2010). The general 
properties of both sources are given in Table 1 . 

3. DATA REDUCTION AND ANALYSIS 



3.1. FERMI-Large Area Telescope 

The Ferm/-Large Area Telescope (LAT; Atwood et al. 
2009) data used in this work were collected over the last 
four years of Fermi operation, from 2008 August 05 00:00:00 
UT to 2012 July 24 00:00:00 UT. LAT photons of event 
class 2 between 100 MeV to 300 GeV were extracted from 
Fermi Science Support Center 4 within 30° circular region of 
interest centered around the position of the objects. Anal- 
ysis of the data was done with the LAT software package 
ScienceTools v9r27pl along with the use of post- 
launch instrument response functions P7SOURCE_V6 and 
the corresponding Galactic and isotropic diffuse background 
models. Additionally, a cut on the zenith angle (< 100°) was 
applied in order to exclude the contamination from Earth's- 
albedo. The spectral analysis were performed using an un- 
binned maximum likelihood method. The data were fitted 
with a simple power law model. In the analysis, we consider 
the source to be detected if test statistic (TS) > 9. This corre- 
sponds to —3cr detection (Mattox et al. 1996). ForTS < 9, 2a 
upper limits were calculated. Systematic uncertainties in the 
flux were estimated as 10% at 0.1 GeV, 5% at 560 MeV and 
20% at 10 GeV and above 5 . Results of the average analysis 
are given in Table 1 . 

3.2. SWIFT (BAT, XRT, UVOT) 

We have searched for the data from all the three tele- 
scopes onboard Swift (Gehrels et al. 2004), namely the Burst 
Alert Telescope (BAT; Barthelmy et al. 2005) operating in the 
15 — 150 keV band, the X-ray telescope (XRT; Burrows et al. 
2005) operating in the 0.3-10 keV band and the Ultra- Violet 
Optical Telescope (UVOT; Roming et al. 2005) which can ob- 
serve in 6 filters; namely V, B, U, UVW1, UVM2 and UVW2. 
We extracted all publicly available optic al/UV/X-ray data for 
the sources from HEASARC archives. The data were analyzed 
using latest Swift tools included in HE AS OFT v . 6 . 1 2 . to- 
gether with calibration data files (CALDB) updated on 2012 
April 02. 

Both 7-NLSyls could not be detected in the 15-150 keV 
band. We analyzed the XRT data using the XRTDAS software 
package developed at the ASI Science Data Center (ASDC) 
and distributed by HEASARC within the HEASoft package 
(v.6.12.0). Data were cleaned with standard filtering criteria 
and calibrated using the xrtpipeline task with default pa- 
rameters and single-to-quadruple events (grades 1 — 12). We 
did not find any evidence of pile-up. Events for spectral anal- 
ysis were selected within a circle of 50" radius, centered on 
the source position, while the background was selected from 
a nearby source-free region with 100" radius. To generate 
lightcurves we first rebinned the spectrum to have at least 1 
count per bin (using grppha) and then calculated the like- 
lihood using C-statistics (Cash 1979). The data were fitted 
with absorbed power-law model with Galactic absorption val- 
ues taken from Kalberla et al. (2005). For lightcurve analysis, 
we have not corrected the XRT data for galactic absorption. 
For SED analysis, the individual XRT event files were com- 
bined together using the XSELECT package and the average 
spectrum was extracted from summed event files. Depending 
on the total detected counts, we grouped the spectrum so as 
to have at least 5 counts per bin for PKS 1502+036 and 15 
counts per bin for PKS 2004-447. Model fitting was done 

4 http://fermi.gsfc.nasa.gov/ssc/data 

5 http://fermi.gsfc.nasa.gov/ssc/data/analysis/LAT_caveats.html 
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Table 1 

List of the 7-NLSyl galaxies monitored in this work and the results of the analysis of about four years of Fermi LAT data. Column informations are as follows: 
(1) Name; (2) right ascension; (3) declination; (4) redshift; (5) absolute B magnitude; (6) apparent V magnitude; (7) radio spectral index; (8) radio loudness 
parameter; (9) 0.1—300 GeV flux in the units of 10 — 8 ph cm -2 s — 1 ; (10) photon index obtained from Ferair'-LAT data analysis; (11) TS and (12) Nh 



Name 


RA (2000) 


Dec (2000) 


z a 


M s a 


V" 


a R h 


R c 


Fo.l-300GeV 


r 


TS 


AV 




(h m s) 


(d m s) 




(mag) 


(mag) 














(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


PKS 1502+036 


15:05:06.5 


+03:26:31 


0.409 


-22.8 


18.64 


0.41 


3364 


5.15 ±0.42 


2.67 ± 0.06 


419.28 


3.93 


PKS 2004-447 


20:07:55.1 


-44:34:43 


0.240 


-21.6 


19.30 


0.38 


6358 


1.66 ±0.33 


2.58 ±0.12 


67.67 


3.17 



a Veron-Cetty & Veron (2010) 

b Radio spectral index calculated using 6 cm and 20 cm flux densities given by Veron-Cetty & Veron (2010) 
Toschini (2011) 

d Galactic absorption in units of 10 20 cm -2 from Kalberla et al. (2005) 



using XSPEC (Arnaud 1996). Finally, the SED data were cor- 
rected for Galactic absorption and then binned appropriately 
for both sources. 

Swift-XJVCfT has observed these sources in all six filters for 
most of the observations but not always. The UVOT data were 
integrated with uvotimsum task and then analyzed with 
uvot source task. The source region was chosen as a circle 
of 5" radius for optical filters and 10" radius for UV, while the 
1 ' sized background region was extracted from nearby source- 
free regions. The observed magnitudes were de-reddened us- 
ing the extinction laws of Cardelli et al. (1989) and converted 
to flux units using the zero points and conversion factors of 
the Swift-CALDB (Poole et al. 2008). 

4. RESULTS 
4.1. Multi-band temporal behaviour 

The long-term multi-band light curves of the two sources 
PKS 1502+036 and PKS 2004-447 are shown in Figure 1 
and 2 respectively. In the 7-ray band, flux values were ob- 
tained over 30 day binning, whereas in the other bands, each 
flux point corresponds to one observation of Swift. For PKS 
1502+036, only four sets of Swift data were available. As the 
data are sparse, no strong claims of variability could be made. 
However, visual examination show some hints of variability 
in the UVOT bands. In the 7-ray band, the source was mostly 
in the low activity state, however, was found to be in a rela- 
tively bright state between MJD 55592 and 55622. It is also 
found to be variable at 15 GHz when observed from Ovens 
Valley Radio Observatory (OVRO) 6 . 

For PKS 2004—447, the sparseness of the data precluded 
us to provide any estimates of the variability of the source. 
Visual examination of the lightcurves in the UVOT bands 
give hints that the source has shown variability. Source has 
also shown significant variability in the X-ray band. When 
observed first time on MJD 55696, its observed flux was 
(4.42±1.23) x 10 -13 erg cm -2 s _1 , which increased by ~3 
times when observed again on MJD 55821 but diminished to 
(5.36±2.20) x 10~ 13 erg cm~ 2 s~ J when observed on MJD 
56130. In the 7-ray band, the source was detected only in two 
epochs when analyzed using a 30 day binning. Comparing the 
counts at those two epochs we find that the source has varied 
in flux. 

4.2. X-ray spectral fit 

For PKS 1502+036, all Swift/XRT observations were inte- 
grated (global exposure = 19.25 ksec) and then fitted with an 
absorbed power law. We find photon index, T = 1 .70±0.20, 

6 http://www.astro.caltech.edu/ovroblazars/data/data.php 



the normalization at 1 KeV = (4.49±0.67) x 10 -5 ph cm~ 2 
s _1 keV -1 and the integrated flux in the 0.3 — 10 keV energy 
band = (2.90±0.69) x 10" 13 erg cm~ 2 s _1 . 

On fitting an absorbed power law model to the integrated 
data (global exposure = 45.89 ksec) of PKS 2004-447, 
we find V = 1.53±0.09 and the normalization at 1 keV = 
(8.14±0.57) x 10" 5 ph cm" 2 s" 1 keV" 1 . The flux in the 
0.3-10 keV energy band = (6.26±0.73) x 10" 13 erg cm" 2 
s- 1 . 

4.3. 7-ray spectral variability analysis 

In Figure 3, we show the correlation of the 0.1—300 GeV 
flux against the 7-ray spectral slope for the source PKS 
1502+036, when 30 day time binning was used. The photon 
index values were obtained by fitting a power law model to the 
Fermi-hAT data. Clear evidence of spectral softening when 
brightening was found for this source when the data was fit 
using both the unweighted least squares method and weighted 
linear least squares fit that takes into account the errors on 
both the dependent and independent variables (Press et al. 
1992). In the case of PKS 2004-447, we have only two 
detections at the epochs MJD 55232 and MJD 55802. The 
flux values at these two epochs are (7.93 ± 2.74) x 10~ 8 and 
(4.72 ± 1.91) x 10~ 8 ph cm" 2 s" 1 respectively. The pho- 
ton indices at MJD 55232 and MJD 55802 are 3.12 ± 0.38 
and 2.56 ± 0.31. Thus, between the two epochs of detection, 
the 7-ray spectral index shows a softening when brightening 
trend. This behavior of the photon index is contradictory to 
the harder when brighter trend seen in the brightest FSRQ 3C 
454.3 (Abdo et al. 2010a). 

4.4. Spectral Energy Distribution 

To investigate the physical characteristics of these sources, 
we have built optical/UV to 7-rays SEDs using quasi- 
simultaneous data. For PKS 1502+036, we have used the 
data sets between MJD 55800 and 56130, whereas for PKS 
2004-447, the data sets spanning the period MJD 55680 and 
56130 were used. These time periods are shown as verti- 
cal dashed lines in Figure 1 and 2 for PKS 1502+036 and 
PKS 2004—447 respectively. The derived flux values in op- 
tical/UV, X-ray and 7-rays are given in Table 2. Both the 
SEDs were modeled using a simple one zone leptonic emis- 
sion model described in Sahayanathan & Godambe (2012). 
In this model the emission region is assumed to be a spher- 
ical blob moving down the jet at relativistic speed with bulk 
Lorentz factor V. The emission region is filled with a broken 
power law electron distribution with indices p and q before 
and after the break energy, respectively. The particles lose 
their energy through synchrotron radiation in a randomly ori- 
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ented magnetic field (B) and inverse Compton radiation. The 
target photons for the inverse Compton radiation are the syn- 
chrotron photons (SSC) and an external photon field which 
is assumed to be a blackbody for simplicity. Considering the 
magnetic field to be in equipartition with the relativistic parti- 
cle distribution, the main parameters governing the SED can 
be deduced/constrained using the observed information avail- 
able in optical/X-ray and 7-ray energies. The kinetic power 
of the jet can then be calculated by assuming the inertia of 
the jet being provided by cold protons whose number density 
is equal to that of non-thermal electrons. This model does 
not take into account UV/optical emission from the accretion 
disk and X-ray emission from the hot corona. For the model 
fitting we assume a variability time scale of ~1 day, which is 
typically seen in the 7-ray lightcurves of blazars. 

The resultant model fits to the SEDs for PKS 1502+036 
and PKS 2004-447 are shown in Figure 4. The optical-UV 
flux can be reproduced by synchrotron emission and the X- 
ray flux by SSC emission. The explanation of the 7-ray flux 
as an outcome of EC process demands the temperature of the 
external photon field to be 605 K in case of PKS 1502+036 
and 694 K for PKS 2004-447. The spectrum of this external 
photon field black body component with luminosities of 8. 19 
x 10 43 erg s- 1 and 6.26 x 10 43 erg s" 1 for PKS 1502+036 
and PKS 2004—447 respectively is shown in Figure 4. We 
compared the luminosity of this black body component with 
the observed luminosity (which is the sum of thermal and non- 
thermal components) in 2MASS and WISE bands and as ex- 
pected they are lower than the observed luminosities. Also 
the deduced temperatures are consistent with the temperature 
generally observed for the dusty torus (Jaffe et al. 2004). The 
archival data from NED is shown with green color in Figure 
4. The radio observation may be contaminated by the core 
emission and emissions from other extended region and hence 
do not satisfy the model curve. The model fitting parame- 
ters are given in Table 3. We note that many of the parame- 
ters obtained from our modelling of the broad band SED of 
these two sources differ from the paramaters obtained for the 
same sources by Abdo et al. (2009c). This might be partly 
due to the quality of the data sets used for SED modelling. 
Here we use nearly simultaneous data sets compared to the 
non-simultaneous multi-wavelength data used by Abdo et al. 
(2009c). 

5. DISCUSSION 

Aligned relativistic jets are invoked to explain the extreme 
power generated by the blazar class of AGN, with a large frac- 
tion of the power being emitted in the 7-ray band. RL-NLSyl 
galaxies also have many observed characteristics similar to 
blazars, namely, flat/inverted radio spectra, compact core-jet 
structure on VLBI scales (Giroletti et al. 2011; Orienti et al. 
2012), high brightness temperature, superluminal jet compo- 
nents (Doietal. 2006), INOV (Paliya etal. 2013; Liuetal. 
2010) and rapid infra-red flux variability (Jiang et al. 2012). 
However, the detection of variable 7-ray emission in five 
RL-NLSyl galaxies (Abdo et al. 2009b,c; D' Ammando et al. 
2012; Calderone et al. 201 1) clearly shows that NLSy 1 galax- 
ies could also have relativistic jets similar to the blazars hosted 
in elliptical galaxies. If the NLSyl galaxies are hosted in 
spiral galaxies we should revisit the paradigm that jets can 
only be launched in elliptical galaxies (Marscher 2009). The 
blazar like behavior shown by these 7-NLSyl galaxies, pow- 
ered by smaller mass black holes and accreting close to the 
Eddington limit, also, poses questions on the nature of them 



vis-a-vis blazars and their place in the traditional blazar se- 
quence. Thus an understanding of these sources will enable 
one to probe jet physics at a different mass range compared to 
blazars. We have presented here a detailed analysis of two 
such RL-NLSyl galaxies namely PKS 1502+036 and PKS 
2004—447 using multi-wavelength data and modeling. Our 
multi-wavelength lightcurves clearly show the objects to be 
variable; however, owing to the sparseness of the data, cor- 
relations if any on the flux variations between different fre- 
quency bands could not be ascertained. 

From the modeling of the SEDs of both PKS 1502+036 
and PKS 2004—447, using the one zone leptonic model of 
Sahayanathan & Godambe (2012), we find that, for these 
sources, the UV-optical data are well fit with synchrotron 
emission, the X-ray data are fit with SSC and the 7-ray points 
are fit with EC scattering of IR photons from the torus. The 
power of the jet estimated from the fit parameters is larger 
for PKS 1502+036 compared to PKS 2004-447. To see how 
the SEDs of these two 7-NLSyl galaxies match with that of 
blazars, we also built two other SEDs, one for the well stud- 
ied FSRQ 3C 454.3 (2 = 0.859) and the other for the BL Lac 
Mrk421 (z = 0.030). For both 3C 454.3 and Mrk 421, we uti- 
lized data that covered the same period as that used for the two 
7-NLSyl galaxies studied here. The data for both 3C 454.3 
and Mrk 421 were analyzed using the procedures outlined in 
Section 3 except for the XRT analysis of Mrk 421 wherin we 
have used a log-parabola model (Massaro et al. 2004a,b). The 
results of the spectral analysis for 3C 454.3 and Mrk 421 are 
given in Table 4. In Figure 5, we show the SEDs of PKS 
1502+036 and PKS 2004-447 along with the FSRQ 3C 454.3 
and BL Lac Mrk 42 1 . In this Figure too, the archival data are 
from NED and BZCAT 7 . We found that during the period con- 
sidered for spectral analysis, 3C 454.3 was in quiescent state 
while Mrk 421 was in exceptional bright state. 

From Figure 5, it is clear that for 3C 454.3, the synchrotron 
emission peaks in the IR region, whereas the high energy 
EC emission peaks in the MeV— GeV range. However, for 
Mrk 421, the synchrotron emission peaks in the UV— X-ray 
spectral band and the high energy SSC emission peaks in the 
GeV— TeV region. While the optical/UV luminosities of the 
two 7-NLSyl galaxies are comparable to that of Mrk 421, 
their 7-ray luminosities are significantly higher than the flar- 
ing Mrk 421. This clearly demonstrates that SSC alone can- 
not account for the high energy emission from these two 7- 
NLSyl galaxies and that EC is needed. Such EC model is 
frequently invoked to explain the high energy emission from 
FSRQs (Dermeretal. 1992; Sikoraetal. 1994, 2009). In 
the powerful FSRQ, PKS 1510-089, based on the observed 
coincidence of 7-ray flux to the appearance of new knots, 
Marscher et al. (2010) hint that the 7-rays are produced from 
a region well beyond the radius of the BLR, most likely in the 
IR torus. Thus, the high energy emission from these 7-NLSyl 
galaxies clearly resemble that of FSRQs and likely due to IC 
scattering of IR photons from the dusty torus. Another sup- 
porting fact for the above mentioned theory is that nominal 
temperature of molecular torus is about ~ 300 K (Landt et al. 
2010) to - 1200 K (Cleary et al. 2007; Malmrose et al. 201 1) 
with hot emitting region expected to be nearer to the cen- 
tral super massive black holes (SMBHs). The IR torus tem- 
perature which we obtained from SED analysis is about ~ 
600—700 K for both sources, which means that the dissipa- 
tion region may lie outside the BLR. Recently, Zhang et al. 

7 http://www.asdc.asi.it/bzcat/ 
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(2012) too have noted that the observed broad band SEDs of 
four NLSyl galaxies detected by Fermi are in close resem- 
blance to FSRQs than BL Lacs. In the SEDs of these NLSyl 
galaxies (see Figure 5), it is clear that the low energy bump 
is in the IR regime, while the high energy bump is in the 
MeV— GeV band. Also, it is clear from Figure 5, that both 
the low and high energy peaks of these 7-NLSyl galaxies are 
intermediate to that of 3C 454.3 and Mrk 421. Thus, both 
PKS 1502+036 and PKS 2004-447 (and may be other 7- 
NLSyl galaxies) could fit very well into the traditional blazar 
sequence (Fossati et al. 1998; Foschini et al. 2010). 

Recently, Chen & Bai (2011) has reported that the 7- 
NLSy 1 galaxies occupy low luminosity-low frequency region 
in synchrotron peak frequency (jA^-luminosity (L s ) plane. 
This trend can also be very well seen in Figure 5, where L s of 
both PKS 1502+036 and PKS 2004-447 are similar to that of 
Mrk 421, whereas, their corresponding peak frequencies are 
quite different. Both of them have v s in the lower frequency 
region compared to Mrk 421. Again the possible explanation 
for such low luminosity - low frequency behaviour of both 
PKS 1502+036 and PKS 2004-447 could be due to their low 
black hole masses (Chen & Bai 201 1; Ghisellini & Tavecchio 
2008). 

In Figure 6, we have shown the position of the two NLSyl 
galaxies studied here, in the 7-ray spectral index v/s K- 
corrected 7-ray luminosity plane along with the FSRQ 3C 
454.3 and the BL Lac Mrk 421. Also shown in Figure 
6, are the other three known 7-ray emitting NLSyl galax- 
ies namely 1H 0323+342 (z = 0.063), SBS 0846+513 (z = 
0.584) and PMN J0948+0022 (z = 0.584). The K-corrected 
7-ray luminosities were evaluated following Ghisellini et al. 
(2009). Considering only the sources studied here, namely, 
PKS 1502+036 and PKS 2004-447, we find that their 7- 
ray luminosities are intermediate to the FSRQ 3C 454.3 and 
BL Lac Mrk 421 and, they have steep 7-ray spectral index 
similar to FSRQs (Ghisellini et al. 2009). However, consider- 
ing all the Fermi detected NLSyl galaxies till now, we found 
that they have a wide range of 7-ray luminosities overlapping 
the region occupied by BL Lacs and FSRQs (Ghisellini et al. 
2009, 2011), however, the steep 7-ray spectral index dis- 
played by most of them, leads us to argue that in terms of the 
7-ray spectral index behavior, Fermi detected NLSyl galaxies 
are similar to FSRQs. 

Further, as discussed by Costamante (2009), FSRQs are 
more Compton-dominated (Compton-dominance is the ra- 
tio between Compton and synchrotron luminosities, CD = 
Lc/Lg) than BL Lacs and on average show a sequence of 
values from ~ 100 in FSRQs to < 1 in high energy peaked BL 
Lac objects (HBLs). If we calculate CD for our sources, it 
comes ~ 10 for both of them. Finally, X-ray spectra of low- 
energy peaked BL Lac objects (LBL) and FSRQs are believed 
to be dominated by IC emission of low-energy electrons and 
thus a flat spectral index (a x < 1) is expected. As we move 
from LBL to HBL sources, the peak shifts toward higher en- 
ergies and tail of the synchrotron emission starts to dominate 
and a steep spectrum (a x > 1) is expected. If we calculate the 
X-ray spectral indices for our sources (see Table 2 and 4), we 
found that both PKS 1502+036 and PKS 2004-447 are hav- 
ing flat spectral indices (0.88 and 0.53, respectively), which 
are intermediate to 3C 454.3 (0.49) and Mrk 421 (1.31). 

Thus, from an analysis of the presently available multi- 
wavelength data set of PKS 1502+036 and PKS 2004-447, 
we find that many properties of them are intermediate to the 
FSRQ 3C 454.3 and the BL Lac Mrk 421. Further, our anal- 



ysis also suggests that these sources may be low black hole 
mass counterparts of powerful FSRQs. However, we caution 
that as of now there are no unambiguous observational evi- 
dence for these two sources in particular and NLSyl galaxies 
in general, to be hosted by spiral galaxies having low mass 
black holes. However, going by the present widely consid- 
ered notion that NLSyl galaxies are powered by low mass 
black holes in spiral galaxies, the current observational evi- 
dence for the presence of aligned jets in the NLSyl galaxies 
detected by Fermi goes against the prevailing idea of jets be- 
ing launched only by elliptical galaxies. A clear picture on the 
nature of these sources will emerge when more observations 
become available in the future. 

6. CONCLUSIONS 

In this work, we have analyzed the two intermediate red- 
shifted 7-NLSy 1 galaxies using all available data covering IR, 
optical, UV, X and 7-ray regimes to understand their physical 
characteristics. We also compare the results obtained, with 
two other well studied sources namely the FSRQ 3C 454.3 
and the BL Lac Mrk 421. We arrive at the following conclu- 
sions. 

1. The broad-band SEDs of PKS 1502+036 and PKS 
2004—447 resemble more, although less powerful, to 
the SEDs of the FSRQ class of AGN. A possible in- 
terpretation is that they could be the low black hole 
mass counterparts of powerful FSRQs. However, the 
potential connection between NLSyl galaxies and FS- 
RQs needs further detailed investigation. 

2. The Compton dominance and X-ray spectral index of 
PKS 1502+036 and PKS 2004-447 are between the 
values found for 3C 454.3 and Mrk 421. In the 7-ray 
luminosity — spectral index plane, these sources have 
luminosities intermediate to that of 3C 454.3 and Mrk 
42 1 , however, they have steep 7-ray spectral index sim- 
ilar to the FSRQ 3C 454.3. On inclusion of the other 
three 7-ray emitting NLSyl galaxies (see Figure 6), we 
found that 7-NLSy 1 galaxies may have a wider range of 
7-ray luminosities between FSRQs and BL Lacs, how- 
ever, they have 7-ray spectral index more like FSRQs. 
The derived spectral properties of the two sources stud- 
ied here are, in general, similar to blazars and inter- 
mediate between the FSRQ 3C 454.3 and BL Lac Mrk 
421. They thus, could fit well into the traditional blazar 
sequence. 

3. In the 7-ray band, spectral softening when brightening 
trend is clearly seen for PKS 1502+036. Similar trend 
is also seen in PKS 2004—447 between the available 
two epochs of observations. This is against the harder 
when brighter trend seen in the FSRQ 3C 454.3. 



This research has made use of the data obtained from the 
High Energy Astrophysics Science Archive Research Center 
(HEASARC) provided by the NASA's Goddard Space Flight 
Center. This research has made use of the NASA/IPAC Extra- 
galactic Database (NED) which is operated by the Jet Propul- 
sion Laboratory, California Institute of Technology, under 
contract with the National Aeronautics and Space Adminis- 
tration. This research has made use of the XRT Data Anal- 
ysis Software (XRTDAS) developed under the responsibility 
of the ASI Science Data Center (ASDC), Italy. 



6 



V S Paliya et al. 



Table 2 

Summary of SED analysis 



Fenni-hPiT 


Source 


Time covered 11 


Flux" 


Photou Index L 


Test Statistic" 








PKS 1502+036 


55800-56130 


2.1 ±0.31 


2.93±0.19 


87.17 








PKS 2004-447 


55680-56130 


0.71±0.19 


2.93±0.30 


19.17 








Swift-XBJ 


Source 


Time covered 11 


Exp. L 


r 1 


Flux s 


Normalization 11 


Stat. 1 




PKS 1502+036 


56042-56103 


14.58 


1.88±0.22 


2.88±0.72 


5.28±0.82 


14.26/21 




PKS 2004-447 


55696-56130 


45.89 


1.53±0.09 


6.26±0.73 


8.14±0.57 


36.69/37 




Swift-UVUI 


Source 


Time covered' 1 


v j 


b> 


u 1 


uvwl J 


uvm2' 


uvw2' 


PKS 1502+036 


56042-56130 


6.54±1.26 


5.08±0.82 


4.56±0.48 


4.53 ±0.40 


3.61±0.51 


5.47±0.32 


PKS 2004-447 


55696-56130 


6.59±0.52 


4.84±0.34 


3.35±0.23 


1.62±0.02 


1.18±0.12 


1.21±0.02 



"Time covered for analysis, in MJD. 

b Integrated 7-ray flux in 0.2—300 GeV energy range in the units of 10 -8 ph cm -2 s — x . 
c Photon index calculated from 7-ray analysis. 
d Significance of detection using likelihood analysis. 
c Net exposure in kiloseconds. 
'Photon index of power-law model. 

s Observed flux in the units of 10~ 13 erg cm -2 s _1 , in 0.3—10 keV energy band. 
h Normalization at 1 keV in 10~ 5 ph cm -2 s _1 keV -1 . 
'Statistical parameters:x 2 /dof. 

^Average flux in Swift V, B, U, Wl, M2 and W2 bands respectively, in units of 10 — 13 erg cm -2 s 



Table 3 

Summary of model parameters for the SED of PKS 1502+036 and PKS 2004-447 



Parameter 


Symbol 


PKS 1502+036 


PKS 2004-447 


Redshift 


z 


0.409 


0.240 


Particle spectral index (low energy) 


P 


2.22 


2.10 


Particle spectral index (high energy) 


q 


4.5 


4.00 


Magnetic field (equipartition) 


B 


0.60 G 


0.72 G 


Bulk Lorentz factor 


r 


12 


9 


Comoving emission region size 


R' 


2.86 X 10 16 cm 


2.32 X 10 16 cm 


IR dust temperature in AGN frame 


r 


605 K 


694 K 


IR torus energy density in AGN frame 


u* 


1.02 X 10 -3 erg cm -3 


1.75 X 10~ 3 erg cm -3 


IR black body luminosity in AGN frame 


L* 


8.19 X 10 43 ergs^ 1 


6.26 X 10 43 ergs" 1 


Particle energy density 


U e 


1.43 X 10~ 2 erg cm" 3 


2.06 x 10 -2 erg cm" 3 


Minimum electron Lorentz factor 


1 min 


30 


60 


Break electron Lorentz factor 


7'b 


1870 


1176 


Maximum electron Lorentz factor 


'y max 


1.8 x 10 4 


7.5 x 10 3 


Jet power 


Pi 


3.1 X 10 45 erg s" 1 


8.5 X 10 44 ergs" 1 



Table 4 

Summary of the comparison source analysis. 



Fermi-LAT 


Source 


Time covered 


Flux 


Photon index 


Test Statistic 












3C 454.3 


55680-56130 


4.9±0.29 


2.61±0.06 


742.73 












Mrk421 


55680-56130 


11.4±0.27 


1.79±0.02 


10971.46 












Swift-XRT 


Source 


Time covered 


Exp. 


Nff 


a 




r 


Flux 11 


Normalization 


Stat. 


3C 454.3 


55729-55924 


17.41 


6.58 






1.49±0.05 


(6.54±0.38) x 10~ 12 


(8.43±0.35) x 10" 4 


86.77/74 


Mrk421 


55621-56078 


112.4 


1.91 


2.308±0.001 


0.437±0.004 




(3.42±0.01) x 10" 10 


(8.96±0.01) x 10~ 2 


1177.25/675 


Swift-UVUT 


Source 


Time covered 


v 4 


b 4 


u 4 


uvwl 4 


uvm2 


uvw2 4 






3C 454.3 


55729-55924 


7.13±0.13 


6.94±0.10 


7.46±0.13 


7.63±0.16 


9.28±0.16 


6.59±0.13 






Mrk421 


55621-56078 


132.98±4.34 




156.31±8.24 


167.59±1.04 


183.96±8.48 


169.82±9.91 







Note. — See Table 2 for description of remaining parameters. 



a Log-parabola model parameters used in XSPEC 

b 0.3-10 keV flux, in the units of erg cm~ 2 s _1 

c normalization having units as ph cm~ 2 s~ 1 keV~ 1 

A Swift V, B, U, Wl, M2 and W2 band fluxes in the units of 10 -12 erg cm -2 s" 1 
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Figure 1. Multi-band lightcurves of PKS 1502+036. The OVRO data are in Jansky. UVOT data are in magnitude. XRT fluxes are in units of 10~ 13 erg cm -2 
s — 1 , while LAT 7-ray data has the unit of 10 — 7 ph cm -2 s — 1 . Here, 95% upper limits are shown by vertical arrows. Dashed vertical line shows the time interval 
covered for SED modeling. 
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Figure 2. Multi-band lightcurves of PKS 2004—447. UVOT data are in magnitude. XRT fluxes are in units of 10 12 erg cm 2 s 1 , while LAT 7-ray data has 
the unit of 10 — 7 ph cm -2 s _1 . Here, 95% upper limits are shown by vertical arrows. Vertical dashed lines indicate the period for which SED has been made. 
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Figure 3. 7-ray photon flux v/s photon-index for the source PKS 1502+036. 
Here, dashed line (blue color) shows the unweighted least squares fit and the 
solid line (red color) shows the weighted least squares fit. 

Facilities: Fermi, Swift. 
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Figure 4. SEDs of 7-NLSyl galaxies PKS 1502+036 and PKS 2004-447. Flux due to Fermi-LAT (filled squares), Swift XRT (filled diamonds) and UVOT 
(filled circles) are indicated with red symbols. Vertical arrow shows the upper limit of 7-ray flux. Archival data are shown in green color: 2MASS (upward 
triangles Skrutskie et al. 2006), WISE (open squares; Wright et al. 2010), GALEX (plus sign; Martin et al. 2005) and Radio (filled star; NED). Dotted black line 
shows synchrotron emission. The SSC (pink color) and EC (green color) processes are shown by dashed and dotted-dashed line, respectively. Blue continous 
line is the sum of all radiative components. Around 10 11 Hz, synchrotron self absorption is clearly visible. 
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Figure 5. SEDs of 7-NLSyl galaxies PKS 1502+036 (left panel) and PKS 2004-447 (right panel) along with that of 3C 454.3 and Mrk 421. Spectral data 
points of 7-NLSyl galaxies are shown with blue colors. Data points of 3C 454.3 (upper side of the plot) and Mrk 421 (lower side of the plot) are shown with red 
and green colors respectively. Vertical arrows show the upper limit of Fermi-LAT 7-ray flux. Fermi-hKT data is shown with filled squares, while Swift XRT and 
UVOT data are shown with filled diamonds and filled circles respectively. Archival radio data is shown with filled stars, whereas IR, UV and hard X-ray data are 
shown with filled upward triangles, plus sign and hourglass respectively. Filled downward triangles show the MAGIC (Anderhub et al. 2009) upper limit data of 
3C 454.3 and HAGAR (Shukla et al. 2012) data of Mrk 421. 
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Figure 6. The spectral index (ct 7 ) against the Luminosity (L 7 ) in the energy 
band 0.1-10 GeV. The horizontal grey line marks a 7 = 1.2. 
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